Can Nonlinear Hydromagnetic Waves Support a Self-Gravitating Cloud?

Using self-consistent magnetohydrodynamic (MHD) simulations, we explore the hypothesis that nonlinear MHD waves dominate the internal dynamics of galactic molecular clouds. We employ an isothermal equation of state and allow for self-gravity. We adopt ``slab-symmetry,'' which permits motions $\bf v_\perp$ and fields $\bf B_\perp$ perpendicular to the mean field, but permits gradients only parallel to the mean field. The Alfv\'en speed $v_A$ exceeds the sound speed $c_s$ by a factor $3-30$. We simulate the free decay of a spectrum of Alfv\'en waves, with and without self-gravity. We also perform simulations with and without self-gravity that include small-scale stochastic forcing. Our major results are as follows: (1) We confirm that fluctuating transverse fields inhibit the mean-field collapse of clouds when the energy in Alfv\'en- like disturbances remains comparable to the cloud's gravitational binding energy. (2) We characterize the turbulent energy spectrum and density structure in magnetically-dominated clouds. The spectra evolve to approximately $v_{\perp,\,k}^2\approx B_{\perp,\,k}^2/4\pi\rho\propto k^{-s}$ with $s\sim 2$, i.e. approximately consistent with a ``linewidth-size'' relation $\sigma_v(R) \propto R^{1/2}$. The simulations show large density contrasts, with high density regions confined in part by the fluctuating magnetic fields. (3) We evaluate the input power required to offset dissipation through shocks, as a function of $c_s/v_A$, the velocity dispersion $\sigma_v$, and the scale $\lambda$ of the forcing. In equilibrium, the volume dissipation rate is $5.5(c_s/v_a)^{1/2} (\lambda/L)^{-1/2}\times \rho \sigma_v^3/L$, for a cloud of linear size $L$ and density $\rho$. (4) Somewhat speculatively, we apply our results to a ``typical'' molecular cloud. The mechanical power input requiredComment: Accepted for publication in Ap.J. 47 pages, 13 postscript figures. Report also available at http://cfa-www.harvard.edu/~gammie/MHD.p

Linear and non-linear theory of a parametric instability of hydrodynamic warps in Keplerian discs

We consider the stability of warping modes in Keplerian discs. We find them to be parametrically unstable using two lines of attack, one based on three-mode couplings and the other on Floquet theory. We confirm the existence of the instability, and investigate its nonlinear development in three dimensions, via numerical experiment. The most rapidly growing non-axisymmetric disturbances are the most nearly axisymmetric (low m) ones. Finally, we offer a simple, somewhat speculative model for the interaction of the parametric instability with the warp. We apply this model to the masing disc in NGC 4258 and show that, provided the warp is not forced too strongly, parametric instability can fix the amplitude of the warp.Comment: 14 pages, 6 figures, revised version with appendix added, to be published in MNRA

Three-Dimensional Simulations of Magnetized Thin Accretion Disks around Black Holes: Stress in the Plunging Region

We describe three-dimensional general relativistic magnetohydrodynamic simulations of a geometrically thin accretion disk around a non-spinning black hole. The disk has a thickness $h/r\sim0.05-0.1$ over the radial range $(2-20)GM/c^2$. In steady state, the specific angular momentum profile of the inflowing magnetized gas deviates by less than 2% from that of the standard thin disk model of Novikov & Thorne (1973). Also, the magnetic torque at the radius of the innermost stable circular orbit (ISCO) is only $\sim2$ of the inward flux of angular momentum at this radius. Both results indicate that magnetic coupling across the ISCO is relatively unimportant for geometrically thin disks.Comment: 4 pages, 4 figures, ApJL accepte